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STATISTICAL  ASPECTS  OF  ICE  GOUGING 
ON  THE  ALASKAN  SHELF  OF  THE 
BEAUFORT  SEA 

W.F.  Weeks,  P.W.  Barnes,  D.M.  Rearic  and  E.  Reimnitz 


INTRODUCTION  was  heavily  marked  by  long  linear  depressions, 

which  we  will  refer  to  as  gouges,  produced  by  the 
A  survey  of  the  bathymetry  of  the  Beaufort  Sea  ploughing  action  of  ice.  The  depths  and  widths  of 

shows  that  large  areas  of  this  marginal  sea  of  the  gouge  incisions  in  the  sea  floor  reached  several 

Arctic  Ocean  have  water  depths  of  less  than  60  meters  and  several  tens  of  meters  respectively, 

meters.  In  this  region  ungrounded  pressure  ridge  with  gouges  occurring  as  both  individual  isolated 

keels  may  protrude  downward  for  at  least  50  m,  events  and  as  multiple  events,  presumably  pro- 

and  ice  floes  containing  such  keels  drift  in  a  gener-  duced  by  projections  on  the  same  pressure  ridge 

al  pattern  from  east  to  west  along  the  Beaufort  keel  gouging  the  sea  floor  as  part  of  the  same  ice 

coast.  Therefore  it  is  reasonable  to  presume  that  motion  (Kovacs  and  Mellor  1974,  Remind/  and 

such  sea  ice  masses  could  interact  with  the  sea  Barnes  1974). 

floor.  Indeed,  ice-related  disturbances  of  the  sea  In  this  paper  we  discuss  some  random-appear 

floor  have  been  inferred  for  some  decades  from  ing  aspects  of  ice-produced  gouges  that  occur 

observations  of  sea  floor  sediments  entrained  in  along  a  190-km  stretch  ot  the  coast  ol  the  Beau 

obviously  grounded  ice  masses  (Kindle  1924).  At  fort  Sea  between  Smith  Bay  and  C  amden  Bay  (re 

the  time,  such  processes  were  largely  of  academic  fer  to  Fig.  I ).  We  also  include  a  brief  discussion  of 

interest,  and  there  was  little  motivation  to  explore  the  statistical  concepts  and  techniques  that  are 

mem  systematically.  utilized.  As  much  of  the  study  area  is  part  ot  the 

With  the  discovery  of  oil  and  gas  along  the  mar-  1979,  1982,  and  1984  lease  sales  offered  by  the 

gins  of  the  Beaufort  Sea  at  Prudhoe  Bay  and  off  State  of  Alaska  and  the  Federal  Government,  we 

the  Mackenzie  Delta,  processes  modifying  the  believe  that  the  results  reported  here  are  o!  mime 

floor  of  the  Beaufort  Sea  became  of  interest  due  to  diate  interest  to  the  engineering  community  in 

their  possible  effect  on  offshore  design  and  opera-  solved  in  offshore  design  for  the  Beaufort  Sea  as 

tions.  Early  side-scan  sonar  and  precision  fathom-  well  as  of  long-term  interest  to  the  scientific  coni 

etry  data  coupled  with  diving  observations  (Shear-  munity  interested  in  near-shore  processes  m  shal 

er  et  al.  1971,  Kovacs  1972,  Pelletier  and  Shearer  low,  ice-covered  seas.  The  paper  concludes  bv  dis 

1972)  showed  clearly  that  much  of  the  sea  floor  cussing  some  of  these  potential  applications 


f  inure  1.  Pur!  of  the  Alaskan  coastline  of  the  Beau  fort  Sea. 
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ENVIRONMENTAL  SETTING 

Because  of  their  importance  to  offshore  design 
in  arctic  areas,  ice-produced  gouges  have  been  the 
subject  of  a  number  of  investigations,  especially 
since  they  were  recognized  as  a  recurring  sea-floor 
feature  in  the  shallow  portions  of  ice-covered  seas. 
Rather  than  review  this  literature  here,  we  will 
simply  mention  publications  of  general  interest 
that  provide  more  exhaustive  reference  lists.  Re¬ 
views  of  early  work  can  be  found  in  Kovaes  (1972) 
and  Kovaes  and  Mellor  (1974).  Early  studies  off 
the  Mackenzie  Delta  are  described  by  Shearer  et 
al.  (1971),  Pelletier  and  Shearer  (1972),  and  by 
Kovaes  and  Mellor  (1974).  Early  work  off  the 
Alaskan  coast  is  reported  by  Skinner  (1971), 
Reimnitz  et  al.  (1972,  1973),  Barnes  and  Reimnitz 
(1974),  and  Reimnitz  and  Barnes  (1974),  More  re¬ 
cent  work  is  discussed  by  Shearer  and  Blasco 
(1975),  Hnatiuk  and  Brown  (1977),  Reimnitz  et  al. 
(1977a, b;  1978),  Barnes  et  al.  (1978),  and  Barnes 
and  Reimnitz  (1979).  These  studies  provide  a  de¬ 
scription  of  the  nature  of  the  gouges,  the  charac¬ 
teristics  of  the  ice  involved  in  the  gouging  process, 
the  general  distribution  of  gouging  along  the 
coast,  and,  to  some  extent,  the  forces  involved  in 


the  process  and  the  rates  of  gouge  recurrence.  In 
most  of  the  studies  little  attention  was  paid  to  how 
the  observed  gouge  parameters  varied  or  to 
methods  for  estimating  infrequent  gouging  events, 
such  as  the  formation  of  deeper  gouges.  Excep¬ 
tions  to  this  are  the  papers  by  1  evvis  ( 1977a. b)  and 
Wahlgren  (1979a, b).  which  examine  the  statistical 
aspects  of  the  gouges  located  in  the  general  urea  of 
the  Mackenzie  Delta. 

Present  evidence  suggests  that  the  Beaufort  Vu 
shelf  has  been  relatively  stable  during  the  last 
10, (XX)  years  (i.e.  major  tectonic  or  glacio-isostatic 
adjustments  have  not  taken  place  (Hopkins 
1967]).  As  sea  level  has  risen  approximately  35  in 
in  this  time  period,  the  entire  sea  floor  of  the  pres¬ 
ent  study  area  was  land  in  the  geologically  recent 
past.  The  gentle  slope  of  the  present  land  surface 
continues  northward  to  a  water  depth  of  60  to  ~0 
m,  where  the  shelf  break  occurs  (Baines  and 
Reimnitz  1974).  f  igure  2  gives  generalized  bathy¬ 
metry  for  the  study  area.  The  broad,  gently  slop 
ing  shelf  is  quite  evident.  When  the  sea-floor  top¬ 
ography  in  the  study  area  is  examined  in  more  de¬ 
tail  it  is  found  to  be  very  complex  (see  map  given 
as  Appendix  A;  Barnes  and  McDowell  1978).  The 
most  notable  features  arc  a  number  of  submerged 
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f  igure  2.  Generalized  bathymetric  chart  oj  the  study  area. 


shoals  and  bottom  irregularities  that  have  been  re¬ 
lated  to  ice  zonation  (Reimnitz  et  al.  1977b).  On 
the  scale  of  the  gouging  it  is  even  more  complex 
(not  shown).  Holocene  sediments  (chiefly  poorly 
sorted  silty  clays  and  sandy  muds )  exhibit  maxi¬ 
mum  thicknesses  of  5  to  10  m  over  the  inner  shelf. 
The  seabed  of  the  region  is  characterized  by  ex 
l rente  diversity  and  variability  of  sediment  types, 
seabed  character,  and  sedimentary  structures. 
Sedimentary  structures  are  dominated  by  wave- 
arid  current-related  processes  inshore  of  10  m,  by 
ice-,  wave-,  and  current-related  processes  between 
10  and  20  m,  and  by  primarily  ice-related  pro¬ 
cesses  out  to  water  depths  of  50  m  or  more,  w  here 
water-related  depositional  processes  again  domi¬ 
nate.  Noteworthy  is  the  nearly  ubiquitous  occur¬ 
rence  of  stiff  silty  clays  in  outcrops  on  the  inner 
shelf. 

The  oceanographic  regime  of  the  region  has 
been  little  studied.  The  near-shore  circulation  ap¬ 
pears  to  be  strongly  wind-driven  during  the  sum¬ 
mer,  with  flushing  rates  and  currents  closely  re¬ 
lated  to  local  winds.  The  most  striking  oceano¬ 
graphic  events  are  waves,  currents,  and  surges  re¬ 
sulting  from  late  summer  storms.  Local  sea-level 
rises  of  3  m  coupled  with  3-m  waves  have  been  ob¬ 
served.  Limited  data  during  the  summer  suggest  a 
general  westward  water  motion  produced  by  the 
prevailing  easterlies,  but  wind-driven  reversals  are 
not  rare.  During  the  winter  the  dominant  currents 
on  the  inner  shelf  are  believed  to  be  the  results  of 


thermohaline  drainage  out  of  the  t  lore 
regime  of  dense,  cold,  salt-rich  water  prv  '<> 
the  formation  of  sea  ice  (Mathews  1981, 

The  ice  regime  of  the  region  shows  great  s nanc¬ 
es  with  season  and  distance  from  shore.  During 
the  summer,  ice  conditions  are  extremely  variable. 
Much  of  the  study  area  is  commonly  ice-free,  with 
the  southern  edge  of  the  multiyear  pack  ice  occur¬ 
ring  between  10  and  100  km  offshore.  New  ice 
starts  to  form  in  October;  during  the  early  stages 
of  its  formation  ice  movement  velocities  nearshore 
are  similar  to  velocities  offshore  (5  km  day  on  the 
average  with  highs  of  35  km  day  during  storms 
[Thomas  and  Pritchard  1979]).  As  the  new  ice 
thickens,  v  elocities  decrease  at  nearshore  locations 
until  the  ice  becomes  truly  fast,  experiencing  typi¬ 
cal  motions  of  only  a  few  tens  of  meters  over  the 
remainder  of  the  winter  (on  occasion  motions  ex¬ 
ceeding  1000  ni  have  been  observed  (Kovacs 
1979]).  At  offshore  locations,  motions  also  de¬ 
crease  somewhat  but  average  movements  still  re¬ 
main  significant  (1  to  2  km  day").  At  times  the 
whole  ice  pack  may  be  nearly  motionless  for  sev¬ 
eral  days.  Numerous  pressure  ridges  form  in  the 
moving  ice,  and  in  shallower  areas  many  of  these 
ridges  become  grounded.  Areas  of  particularly 
heavy  groundings  occur  off  the  barrier  islands  in 
water  depths  of  roughly  15  to  20  m.  In  areas  such 
as  Harrison  Bay  that  are  not  protected  by  barrier 
islands,  large  grounded  ridges  occur  in  shallower 
waters  (roughly  10-m  depth).  Once  the  grounded 


3 


tidge  or  Mannikin  /one  dm  clop-,  ilic  ice  -hoic 
ward  of  tin-  Icalnrc  remain-  [datively  motioiile— 
u lit 1 1  spring 

During  -print:,  which  on  the  eoa-t  ol  the  Beau 
for l  Sea  occur-  in  lime,  melting  allow-  lormeiK 
bottom-fa-t  ice  near  the  -Imre  to  float  I  In-  allow- 
t he  near-hore  ice  to  become  mobile  once  again. 
Main  example-  of  ice  pile-up  and  o\ei  tide  on 
bcache-  are  known  to  tune  occurred  dining  tln- 
period  (Kinac-  and  Sodln  1480),  Within  the  con 
-traint-  pre-ented  by  the  coa-t  and  In  grounded 
ridge-  and  rubble  field-,  the  near-hore  ice  remain- 
mobile  throughout  the  -unnner  mile--  it  di-np- 
pear-  b>  melting  or  i-  blown  out  to  -ea.  Howe'er, 
the  ntawne  area-  of  grounded  ridge-  and  nibble 
often  remain  grounded  throughout  mo-i  or  even 
all  of  the  summer  (Kovac-  14'h,  Bat  tle- and  Reim- 
nit/  14"4).  A— ociated  with  the-e  grounded  ice  fea¬ 
ture-  at  the  18-  to  20- m  water  depth  i-  a  break  in 
the  -eabed  -lope  and  change-  in  gouge  character 
and  -ediment  texture  (Reimmt/  and  Battle-  I9'4) 
Additional  information  on  the  oceanography  and 
-edimcntologv  of  the  -tudy  area  can  be  found  m 
Al’O  ( lh'8). 


DATA  (Dl  l  fdlON 
AND  TKRMINOI  ()(i\ 

Seven  year-  of  data  obtained  between  14"2  and 
lh’4  (excluding  1474)  were  u-ed  in  the  pre-ent 


-tudv.  with  a  total  -ample  lia-khnc  length  ot  ap 
pi oximatelv  I '(to  km  Data  wete  collected  liom 
the  te-eatch  ve-el-  /  "<"i  and  Katiuk  u-mg  a  -ide 
-an  -on. ii  and  a  pieci-mn  lathometei  (2<K)  kll/l 
Both  -v-tem-  wete  eap.ible  ol  le-oivmg  bottom 
leliel  ol  le—  than  Id  cm,  I  he  -idc-caii  iccotd- 
coveted  the  -ea  tlooi  beneath  the  -hip  in  2txi  m  in 
2'0-in  -wath-,  depending  on  --ale  -election.  I  ig 
me-  ’  and  4  ate  a  icpiv-cni.uivc  -onogiaph  a n -I  a 
lathogram  te-pectiv  cly .  I  he  tia-k-  weie  -paced  to 
provide  (airly  evenly  di-(nbuied  -atnpling  along 
the  coa-t  between  Smith  Bay  and  Camden  Bay 
Data  were  obtained  both  m-ide  and  -eaward  ot  the 
barrier  i-land-  tv'  the  2.x -in  i-obath.  I  igute  5  -how- 
the  location-  of  tile  difterent  -ampling  line-.  I  he 
t rackline  navigatii'ii  wa-  plotted  ill  l-km  -eg 
mem-.  I  he  -onogiaph-  and  lathogiam-  were  al-o 
divided  inti'  l  km  -egment-  and  conelated  ditectlv 
with  the  navigation.  A-  -ome  a-pect-  v'l  the  data 
interpretation  are  subjective,  all  the  counting  and 
mea-uring  vva-  performed  by  one  individual  (D. 
Keane)  tv>  minuni/e  vatiation- 

A  complete  ice  gouge  data  reev'td  -licet  -bowing 
all  mea-iiicment-  i-  given  by  Reanc  et  al.  <14M| 
\  de-cnptiv'ii  v'l  the  general  lechm-jue-  that  wete 
u-ed  in  analy/in;’  the  -onograph-  and  fatlmgiam- 
i-  given  in  Ba me-  et  al.  (I4~8I.  but  a  lew  impoi- 
taut  point  -  affecting  the  parameter-  u-ed  in  tlu 
pre-ent  -tudv  -hould  be  mentioned: 

•I  uTuev  wilier  depth  (;).  I  hi-  wa-  determined 
by  averaging  the  water  depth-  ob-erved  at  the  -tan 
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/ -'mure  1.  Sonoiiraph  of  leceouiied  sea  floor.  M  ater  depth  is  20  m 
Record  taken  20  km  \/  of  (  ape  Hulked 
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and  at  the  end  of  eaeti  I -km  section:  as  ■  changes 
are  iix.mll>  gradual  and  reasonably  smooth, 
should  he  a  reasonable  approximation  to  a  spatial- 
h  integrated  value. 

Dominant  eouee  orientation  (ft).  Templates 
were  used  to  remove  horizontal  exaggeration  tront 
the  sonographs  and  to  obtain  all  measurements  ot 
the  estimated  dominant  orientation  to  within  5 
true.  It  should  be  noted  that  the  gouge  orienta¬ 
tions  within  each  line  segment  are  variable  (see 
I  ig.  .'). 

Spatial  gouge  tn\pu’ne)  (\  t.  In  determining 
the  number  ot  gt'uges  per  kilometer  ot  sampled 
traek  (  V).  every  feature  on  the  tathograph  pre¬ 
sumed  to  result  from  lee  contact  with  the  bottom 
was  counted,  including  individual  gouges  pro¬ 
duced  b>  ditterent  segments  of  what  was  probably 
the  same  pressure  ridge  keel  (our  interest  is  in  the 
number  of  gouges  in  the  bottom,  not  in  the  num¬ 
ber  of  ice  events);  these  V  values  were  (hen  cor¬ 
rected  in  order  to  estimate  .V  ,  the  expected  num¬ 
ber  of  gouges  tltat  would  have  been  seen  on  a  1-km 
sampling  line  il  the  ship's  track  were  oriented  nor¬ 
mal  to  the  dominant  gouge  trend.  This  correction 
was  made  by  tis.ng  A  V  sum.  where  a  is  an 
acute  angle  between  the  slop's  track  and  the  gouge 
orientation.  As  most  gouges  are  oriented  parallel 
to  the  coast  and  the  majority  ot  the  sampling  lines 
weie  roughly  normal  to  the  coast,  these  correc¬ 
tions  were  usually  small.  Gouges  with  depths  of 
le"  than  0.2  m  were  not  counted,  as  it  was  often 
difficult  to  identify  positively  all  of  these  small 
gouges  on  the  fathogram.  Actually,  in  the  original 
data  tabulations  (Rearic  et  al.  1981)  a  value  was 
given  for  the  number  of  gouges  in  the  0  to  ().2-m 
range  that  could  be  distinguished  on  the  sonar 
record.  Although  this  value  can  be  useful,  it 
should  not  be  combined  with  the  data  on  gouges 
deeper  than  0.2  m,  as  it  includes  a  number  of 
gouges  that  do  not  cross  beneath  (he  ship's  track 


I  m ure  4  luthogram  o'  tee-gouged  \ea  Hoot 
Oarer  depin  is  1  ft  m  Reeurd  taken  *V  km  >/  <>t  (  ape 
/  la  Ik  en. 


lie.  that  would  not  appear  on  the  lathometer 
track). 

Gouge  depth  id)  The  depth  was  measured  (on 
the  fathometer  track)  vertically  from  the  level  ot 
the  (presumably  undisturbed)  adjacent  sea  floor  to 
the  lowest  point  in  the  gouge  (see  Fig.  6).  Values 
were  grouped  in  20-cm  class  intervals:  in  some 
cases,  because  of  factors  such  as  ground  swell  and 
wind  chop,  it  was  only  possible  to  determine  the 
number  ot  gouges  that  have  depths  greater  than  a 
specified  value,  so  gouges  less  than  0.2  m  deep 
were  not  considered.  Each  individual  gouge  was 
measured.  The  maximum  gouge  depth  (</max)  ob¬ 
served  in  each  kilometer  of  sample  track  was  also 
determined.  It  should  be  noted  that,  because  of  in¬ 
filling  by  sediment,  the  measured  d  values  are  pre¬ 
sumably  less  than  the  d  values  at  the  time  the 
gouge  formed.  It  is  also  important  to  note  that  the 
gouge  depth  is  almost  always  less  than  the  depth  to 
the  base  of  the  ice  mass  that  produced  the  gouge, 
because  of  ploughed  sediment  sliding  back  into 
the  gouge  as  soon  as  the  keel  moves  on. 

Maximum  gouge  width  ("max).  This  measure¬ 
ment  was  taken  between  the  inside  walls  of  the 
gouge  at  the  level  of  the  undisturbed  surrounding 
sea  floor  (see  Fig.  6);  the  maximum  value  in  each 
kilometer  of  sample  track  was  recorded. 

Maximum  lateral  embankment  heigh:  (Amav). 
This  measurement  is  the  maximum  height  (in  each 
kilometer  of  sample  track)  of  the  embankments  of 
sediment  plowed  from  the  gouges,  measured  rela¬ 
tive  to  the  undisturbed  sea  floor  (see  Fig.  6)  and 
occurring  along  the  margins  of  the  gouges. 

it  should  be  noted  that  v  alues  of  d  and  hmax  are 
determined  purely  from  the  linear  fathometer  pro¬ 
files.  It  is,  however,  known  from  other  data  (sono¬ 
grams,  dive  observations,  and  repetitive  track 
lines)  that  both  gouge  depths  and  lateral  embank¬ 
ment  heights  can  vary  considciablv  along  the 
length  of  a  given  gouge. 


I  linin'  6.  Sehematie  drauini;  oj  a  eouiie  \Iiowihk  the  loeiiinws  of  seven;/  niea\iire 
merits  referred  to  in  the  text. 


As  is  dear  from  comparing  our  terms  with  the 
titles  of  papers  in  our  reference  list,  terminology 
for  ice-induced  sea  floor  features  is  far  from  stan¬ 
dardized.  This  should  not  he  a  problem  as  long  as 
individual  authors  clearly  spell  out  their  usage  of 
specific  terms.  There  is  little  we  can  do  here  to  re¬ 
solve  terminology  disputes.  We  would  simply  like 
to  point  out  that  gouge  and  gouging  in  the  present 
study  correspond  to  scour  and  scouring  in  the 
papers  of  Pelletier  and  Shearer  (1972)  and  I  ewis 
( 1 977a, h)  and  to  score  and  scoring  in  the  papers  of 
Kovacs  (1972)  and  Kovaes  and  Mellor  (1974). 

In  the  analysis  the  data  are  commonly  combined 
into  several  different  groups  based  primarily  on 
geographic  location.  A  given  group  is  referred  to 
by  either  a  geoniorphic  characteristic  common  to 
the  group  or  by  the  name  of  a  geographic  feature 
occurring  within  or  near  the  location  of  the  group. 
Specifically,  the  groups  are: 

a.  l  agoons  (lines  2-4,  14,  15;  .1-7,  8,  9,  12,  1.1, 
14;  5  .1.  12;  8-17,  40,  41) 

b.  l  onely  (lines  7-19,  40,  41,  42) 

c.  Harrison  Bay  (lines  2-19,  5-12.  6-22,  7-15; 
note  that  nearshore  lines  6-20,  21,  21,  24, 


and  25  were  not  used,  as  the  sonat  records  in¬ 
dicated  sand  waves  and  other  features  that 
suggested  extensive  movement  of  bottom 
sediment) 

d.  Jones  Islands  (lines  2  15,  1'.  21.  7-11.  66, 
67,  ’’! ;  9-92;  observations  from  north  of  Spy 
Island  to  the  north  of  the  Midway  Islands) 

e.  McClure  Islands  (lines  19.  10,  11.  11.  14; 
5-1.  4;  7-76;  9-44,  61,  65,  66.  78;  observa¬ 
tions  from  Cross  Island  to  Camden  Bay) 

f.  Jones  Islands  and  Past  (a  combination  of  the 
Jones  and  the  McClure  Islands  data  sets;  i.e. 
all  the  data  seaward  ol  the  barrier  islands 
and  east  of  Harrison  Bay) 

g.  Harrison  Bay  and  Past  (a  combination  of  the 
data  sets  from  Harrison  Bay,  Jones  Islands, 
and  McClure  Islands;  i.e.  all  the  data  sea¬ 
ward  of  the  barrier  islands  except  the  four 
tracks  off  I  onely). 

DATA  ANALYSIS 

In  the  following  live  sections  we  will  analv/e 

our  field  results  concerning  a)  gongs  depths,  b) 
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gouge  orientations.  c)  gouge  1 1 c«.j ncnc> .  d)  gouge 
widths.  and  e)  lateral  embankment  heights. 

(•nugc  depths 

I  o  examine  the  distiibution  ot  gouge  depths,  we 
ptepated  hixiogiamx  ol  gouge  depths  tot  dilteienl 
legions.  I  he  genet  al  natute  ol  these  gtaphs  ap¬ 
peared  to  he  a  deeteasine  exponential  xvuh  a  rapid 
tall  oft  in  the  licquencx  ol  oeeimenee  ot  laiget 
gouges.  \  similar  tendenex  has  been  noted  b\  both 
E  exs  is  (ld~'a,b)  and  U  ahlgren  (id^da.h)  lot 
gouges  occurimg  north  ol  the  Maeken/ie  Delta, 
bin  examination  ol  ibis  data  |l  ewix  14"a)  shows 
that  the  numhei  ol  small  gouges  is  significantly 
less  than  would  oeeur  in  an  exponential  model, 
suggesting  that  some  otltei  type  ol  distribution 
might  also  be  a  possibilitx.  As  will  be  seen,  this 
"non-exponential"  decrease  in  (he  number  ot 
small  gouges  is  not  apparent  m  the  data  set  studied 
in  the  present  paper  We  also  note  that  the  as¬ 
sumption  ol  an  exponential  distribution  ot  initial 
gouge  depths  is,  as  a  first  approximation,  reason¬ 
able  m  that  the  depths  of  pressure  ndge  keels 
measured  to  the  north  ol  the  studs  area  In  sub¬ 
marine  sonar  can  also  he  well  described  In  the  use 


>t.  .n'uqt*  f  ,  tV  > 

figure  7.  Semilug  plot  of  the  number  of  gouges 
observed  vs  gouge  depth  (til  for  four  regions  along 
the  Alaskan  coast  id  the  Beaufort  Sea. 


o*  an  exponential  dixtnbulion  iWadhamx  and 
Home  IdKO). 

1  igute  '  shows  a  sennlog  plot  ol  the  numhei  o! 
gouges  with  dilleieut  gouge  depths  lot  lout  tepte 
sentatise  areas  ot  the  studs  tegion:  I)  liotn  the 
lagoons  (41  data  points).  2)  trom  Haitisoii  lias 
(X42  data  poinlx).  4)  Irom  oil  I  onely  \2s(g)  data 
points),  and  4)  iiom  the  profiles  seassaiu  ol  the 
hairier  islands  and  east  ol  Hatii'on  Has  <16,62(1 
data  points).  ( )t her  groupings  ol  the  data  and  data 
from  other  areas  gaxe  similar  plots  I  tie  lout 
cutxcs  are  xxell  sepaiated  because  the  numhei  ol 
gouges  obseised  in  the  lout  legions  aie  quite  dil 
lerent.  I  his  is  the  result  ol  dittoing  lengths  ol 
sampling  line  and  ol  diflcimg  'palial  gouge  lie 
queuctes  associated  ssiili  shlleieuces  in  ssatei 
depth.  II  the  same  seis  ol  data  ate  plotted  as 
relatixe  Ire  queues  (the  propomon  ol  the  total 
number  ol  observations  I  tom  that  region  that  ox 
cm  s  in  each  ol  the  0.2-m  depth  classes),  the  shapes 
id  die  curses  are  identical  but  tlieie  is  considerable 
overlap.  Note  that  all  plots  are  reasonable  linear 
over  the  complete  range  of  lour  decades  <r  values 
vary  from  D.S>4  to  0.48  | r  gises  the  fraction  ol  the 
sanation  in  ilte  number  ol  gouges  observed  ac¬ 
counted  lot  by  the  regression  line;  in  this  case.  44 
to  *4S"'o | ).  Ibis  suggests  that  the  ulih/aiion  of  an 
exponential  distribution  in  the  Mackenzie  studies 
(lewis  I4~~a)  is  leasonable  as  an  initial  approxi¬ 
mation.  However  we  note  that  even  (hough  the 
correlation  coefficients  are  high,  v  tests  for  good¬ 
ness  of  fit  are  commonly  failed,  future  studies  of 
gouging  should  explore  the  possibility  of  either 
finding  a  more  satisfactory  distribution  function 
or  of  better  rationalizing  the  deviations  from  ex- 
ponentiality . 

I  lie  exponential  distribution  o  a  convenient . 
well-studied  distribution  (see.  among  others,  Hen 
lamin  and  Cornell  1"'0.  Millet  and  1  reund  Id") 
If  the  simple  ftequenex  distribution  is  a  negative 
exponential,  then  the  probability  density  function 
(I’DI  I  of  V  will  also  be  of  a  similar  form 

/yd  l  )  kg  1  •  l>. 

(Here  s  represents  the  values  that  the  tandom  sari 
able  A  may  acquire.)  Because  the  integral  of./yfv) 
from  0  to  <*  must  equal  I.  as  it  contains  all  the 
sample  points  with  non-zero  probabilities. 
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Fills  gives  the  following  l’DI 

7x|x>  Xe  M  x  •  0  (I) 

More  (ho  inuMimim  likelihood  osimiaio  ol  i ho  free 
parameter  X  is  simply  the  reciprocal  of  the  sample 
moan  (\): 

X  1  \ 

Flic  probability  that  a  raiulom  variable  will  as¬ 
sume  a  value  in  the  interval  (v  ,  x  )  is  then 


Ha  •  A  •  v | 


/X(.v)r/.x 


xVx. 


The  cumulative  distribution  function  (C  OT)  is. 
in  turn,  found  by  integration 


/  s(V) 


/’I  A  a  |  j  ,/x(M)r/w 


where  v  •  0.  finally,  because  we  are  interested  in 
the  probability  of  occurrence  of  gouges  that  have 
depths  greater  than  or  equal  to  some  specified 
value,  we  are  largely  concerned  with  the  value  of 
the  exceedance  probability  given  by  the  comple¬ 
mentary  distribution  function  (>\(v) 

<»\(.v>  /'[A  ~.xj  I  /X(x)  e  K>.  (2) 

is  a  particularly  simple  function  to  graph  as 
it  is  a  straight  line  on  semilog  paper  and  has  a 
value  of  I  at  x  0.  Therefore  the  simple  relation 
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can  be  used  to  estimate  /»[/>»</)  (the  expected 
number  of  gouges  with  depths  greater  than  or 
equal  to  d.  given  that  ,V  gouges  have  occurred). 
Values  for  X  for  the  four  data  sets  shown  in  Figure 
7  are  given  in  Table  la.  In  determining  X  the  fact 
that  the  0  to  0.2-m  gouge  depth  class  was  excluded 
was  handled  by  letting  d‘  =  (d  c)  where  c  =  0.2 
m.  the  cutoff  value.  (Note  that  in  Figure  9  the 
nominal  d  =  0  location  is,  in  fact,  d  0.2  m.) 


Note  also  that  when  the  number  of  gouges  is 
fiii'/i,  only  gouges  with  depths  equal  to  or  greater 
than  0.2  m  are  counted.  The  use  of  a  cuiofl  has  an 
undesirable  effect  on  the  estimates  ol  the  mean 
gouge  depth  m  that  the  d  value  obtained  depends 
upon  the  cutoff  in  use  (in  fable  I  the  value  du  , 
refers  to  a  mean  gouge  depth  calculated  using  the 
0.2-m  cutoff).  Io  facilitate  comparisons  between 
our  data  set  and  those  of  other  investigators  we 
also  include  d\ alues  in  Table  I  that  are  calculated 
by  first  estimating  the  number  of  gouges  m  the  0 
to  0.2-m  class  interval  by  exponential  extrapola¬ 
tion  and  then  including  this  estimate  in  the  calcu¬ 
lation  of  the  mean.  I  he  use  of  the  resulting  values, 
of  course,  implicitly  assumes  that  the  distribution 
of  gouge  depths  is  exponential.  I  he  values  given  in 
parentheses  in  Tables  la  and  lb  for  theO.I-m  class 
interval  are  the  extrapolated  values. 

It  is,  however,  possible  to  sharpen  up  the  above 
by  noting  that,  at  least  oil  the  Mackenzie  Delta, 
the  nature  of  the  gouge  depth  distribution  is 
known  to  change  with  water  depth  (I  ewis  1977a). 
This  is  hardly  surprising  in  that  the  large  ice 
masses  that  presumably  produce  the  deep  gouges 
observed  in  deeper  water  are  not,  because  ol 
grounding,  available  to  produce  similar  gouges  in 
shallow  water. 

We  will  now  examine  the  effect  of  such  a  v.iria 
non  within  our  study  area.  That  similar  changes 
will  be  found  to  occur  in  the  Beaufort  Sea  can  be 
surmised  from  Figure  7  in  that  the  shallower  areas 
(lagoons,  Harrison  Bay)  show  no  deep  gouges. 
The  X  and  d  values  corresponding  to  various  5-m 
water  depth  classes  in  the  different  regions  are 
given  in  Table  la,  and  the  X  values  are  plotted 
against  water  depth  (c)  in  Figure  8.  There  is  clearly 
a  general  decrease  in  X  with  increasing  c  w  ithin  the 
range  of  the  data  set. 

For  a  discussion  of  the  area  in  general,  we  have 
combined  all  the  data  for  offshore  areas  unpro¬ 
tected  by  barrier  islands  (lonely,  Harrison  Bay, 
Jones  Islands  and  east)  into  one  data  set  (Table 
lb).  Figure  9  gives  three  representative  plots  of 
data  from  this  combined  set  for  three  different 
5-m  water  depth  intervals,  and  shows  the  fitted 
curves  based  on  eq  1.  Figure  10  shows  the  seven  x 
values  for  this  combined  set  plotted  versus  We 
have  chosen  to  fit  the  %  vs  ;  data  with  a  negative 
exponential  (r  0.95)  purely  as  a  matter  ol  con¬ 
venience,  This  curve  should  not  be  extrapolated 
beyond  the  range  of  the  data.  For  instance,  it  is 
known  (I  ewis  1977a)  that  gouges  off  the  Macken 
zie  Delta  do  not  appear  on  the  sea  floor  at  water 
depths  greater  than  80  m  and  show  a  pe.ik  in  the 
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Table  I  (eonf'd).  Summar;  of  gouge  depth  (</)  measurements. 


ReiQt i ve  F r equency 


Figure  8.  Plot  of  X  (nl  )  vs  water  depth  (/)  in 
meters  fur  four  geographic  areas  along  the 
Alaskan  coast  of  the  Beaufort  Sea. 


d.Gouqe  Depth  (ml 


Figure  9.  Relative  frequency  of  occurrence  oj 
gouges  of  differing  depths  based  on  all  data  from 
offshore  areas  unprotected  by  barrier  islands. 


figure  10.  X  values  (m  )  v.s  water  depth  </)  based 
on  the  data  set  from  o  ffshore  areas  unprotected  by 
harrier  islands. 


mean  gouge  density  at  a  water  depth  of  23  m. 
Therefore,  one  might  expect  that  in  the  present 
study  area  values  may  increase  again  at  ;  •  35 
rn. 

Clearly ,  water  depth  is  a  most  important  param- 
etei  in  studies  of  gouging. 

Couge  orientation 

Determining  the  absolute  cartographic  orienta¬ 
tion  of  every  gouge  would  be  very  time-consum¬ 
ing.  To  provide  some  information  on  gouge  orien¬ 
tations  we  have  visually  estimated  the  dominant 
orientation  that  exists  along  each  kilometer  of 
sample  track.  These  orientation  values,  however, 
do  not  provide  Information  on  the  actual  direction 
of  the  ice  movement;  for  instance,  the  direction 
90  indicates  only  that  the  gouge  runs  along  the 
90"  to  270  '  line  (in  an  east-west  direction). 

Figure  1 1  shows  linear  histograms  of  the  proba¬ 
bility  of  the  occurrence  of  different  orientations. 
The  data  are  displayed  between  0  and  180  .  T  his 
proved  to  be  convenient,  as  there  was  a  natural 
break  in  the  observations  at  this  orientation  (i.e. 
very  few  gouges  were  aligned  north-south).  Sum¬ 
mary  statistics  for  these  observations  are  pre¬ 
sented  in  Table  2.  The  mean  given  here  is  the  cir¬ 
cular  mean  as  calculated  for  axial  data;  the  circu¬ 
lar  variance  has  a  value  near  0  if  the  data  arc  tight¬ 
ly  clustered  and  a  value  near  I  if  the  directions  are 
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TabU-  2.  Descriptive  statistics  on  variations  in  the  dominant  gouge  orienta¬ 
tion. 
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(71  to  83 °D  to  the  coast,  although  this  tendency 
is  not  evident  in  the  measurements  made  off  the 
Jones  Islands.  It  is  reasonable  to  expect  a  ground¬ 
ed  floe  to  rotate  and  move  toward  the  coast  (this 
effect  has  been  observed  in  radar  imagery  at  Bar- 
row  by  Shapiro  (pers.  comm.)).  However,  it  is  not 
clear  to  us  why  this  phenomenon  should  be  more 
pronounced  in  shallow  water.  The  mean  gouge 
orientation  in  the  lagoons  is  99I'7\  which  is  similar 
to  the  gouge  orientations  in  water  deeper  than  2(1 
m. 

Three  factors  presumably  control  the  orienta¬ 
tion  of  the  gouges.  The  first  factor  is  the  wind  di¬ 
rection,  which  at  Kaktovik  is  predominantly  in 
two  directions:  from  the  F.NF.-E  (55 I00‘,7') 
35%  of  the  time  and  WSW-W  (235"-280°D  23% 
of  the  time  (the  mean  wind  speed  is  the  same  (6.7 
m  s|  in  both  directions  (APO  I978|).  Therefore 
the  ice  drift,  which  is  roughly  45  to  the  right  of 
the  surface  wind,  would  be  expected  to  be  between 
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lot)  and  145  7.  a  value  range  just  above  that  ob¬ 
served.  Secondly,  because  the  fast  ice  edge  gen¬ 
erally  parallels  the  isobaths,  ice-ice  interactions 
tend  to  force  the  near-shore  ice  to  drift  parallel  to 
the  coast  even  when  the  tree-drift  direction  is  not 
exactly  parallel  with  the  coast. 

Finally,  as  mentioned,  the  higher  resistance 
commonly  encountered  by  grounded  features  on 
the  near-shore  (shallow  water)  side  of  a  Hoe  will 
cause  gouges  to  form  at  angles  less  than  expected 
from  a  free-drift  situation.  The  end  result  is  there¬ 
fore  an  average  gouge  direction  in  the  range  of  80 
to  125  T,  as  observed. 

Gouge  frequency 

We  now  have  a  reasonable  description  of  the 
probability  of  a  gouge  having  different  gouge 
depths,  given  that  a  gouge  has  occurred.  Next  we 
need  to  determine  how  many  gouges  have  occurred 
so  that  we  can  estimate  ,V  in  eq  3.  The  number  of 
gouges  that  is  of  primary  interest  is  the  temporal 
gouge  frequency,  that  is,  the  number  ot  gouges 
that  intersect  a  unit  length  of  line  per  unit  of  time 
(gouges  per  kilometer  per  year).  As  will  be  seen, 
data  leading  to  such  estimates  are  extremely 
sparse.  What  is  available  are  measurements  of  the 
spatial  gouge  frequency  (e.g.  gouges  per  kilom¬ 
eter)  as  seen  at  a  given  location  at  essentially  a 
fixed  instance  in  time.  We  will  now  discuss  these 
two  parameters. 

Spatial  ftoiifU’  frctpu’tuy 

To  study  variations  in  the  spatial  gouge  fre¬ 
quency.  the  number  of  gouges  deeper  than  0.2  m 
per  kilometer  was  determined  for  each  kilometer 
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I  inure  12.  S' umber  of  gouges  per  kilometer  measured  normal  lo  the 
trend  of  the  gouges  fSd  vs  water  depth  It). 


of  sampling  track.  These  values  were  then  con¬ 
verted  to  ,Vi,  the  number  of  gouges  per  kilometer 
that  would  have  been  encountered  if  the  sampling 
track  were  oriented  perpendicular  to  the  trend  of 
the  gouges.  The  values  were  then  separated  into 
five  different  groups  (lagoons.  Lonely,  Harrison 
Bay,  Jones  Islands,  and  McClure  Islands  and  east) 
and  plots  were  made  of  S',  vs  water  depth.  Exam¬ 
ination  of  these  plots  showed  that  lagoons  differed 
front  the  other  four  areas  in  that  gouges  were  rare 
(92®'o  of  the  298  kilometers  sampled  contained  no 
gouges  and  the  largest  S',  value  was  12  gouges 
km).  The  four  other  regions  showed  differences 
but  these  appeared  to  be  largely  caused  by  changes 


in  the  water  depths  sampled  in  the  different  areas. 
Therefore  all  four  regions  were  combined  and 
considered  as  one.  Figure  12  shows  the  .V,  vs  ;  plot 
for  the  combined  data.  A  data  tabulation  is 
presented  in  Table  3.  As  was  the  case  in  the 
lagoons,  in  shallow  water  ,V,  values  of  zero  (.V,, 
values)  are  common  and  ,V,  values  greater  than  50 
are  rare.  In  water  15  to  20  tn  deep,  zero  values  be¬ 
come  less  common  and  larger  V  values  are  en¬ 
countered.  Finally,  as  water  depths  increase  above 
22  nt,  all  samples  show  20  or  more  gouges  per  kil¬ 
ometer. 

These  changes  can  be  shown  (I  ig.  13)  by  taking 
10-m-wide  vertical  slices  through  Figure  12  and 


I  able  3.  Summary  of  the  number  of  gouge*  per  kilometer  deeper  than  0.2  m. 
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displaying  the  results  as  histograms  giving  relative 
frequency  vs  V,  10.  As  can  be  seen,  in  the  lagoons 
there  is  a  rapid  exponential  drop-off  in  frequency 
as  the  V  10  value  increases.  In  shallow  water 
( 10  in)  outside  of  the  barrier  islands  the  trend  is 
similar,  although  null  values  are  not  as  frequent 
(42ro).  At  depths  of  10  to  20  m,  the  null  values 
compose  only  24°'o  of  the  total  sample  and  .V  10 
values  in  excess  of  10  are  not  rare.  In  deeper  water 
the  distribution  from  20  to  30  m  and  30  to  38  m 
had  nearly  identical  means  and  forms  and  were 
therefore  combined.  The  histogram  is  now  more 
clearly  Gaussian  and  shows  only  a  slight  positive 
skew.  Again,  the  nature  of  the  distribution  is 
clearly  a  function  of  the  water  depth.  One  addi¬ 
tional  piece  of  information  should  be  added  here: 
at  one  location  (off  Lonely)  a  study  was  made  of 
the  distribution  of  the  spacings  between  the 
centers  of  gouges  (as  measured  along  the  sampling 
line);  the  distribution  resembled  a  negative  ex¬ 
ponential  (Fig.  14). 


It  would  be  convenient  to  have  one  distribution 
function  that  would  describe  all  the  histograms 
shown  in  Figure  13.  If  possible  this  distribution 
should  have  the  following  characteristics: 

a.  It  should  be  discrete  in  that  we  are  describing 
a  counting  process  (either  a  gouge  is  present 
or  it  is  not) 

b.  It  should  be  capable  of  dealing  with  a  posi 
live  frequency  of  zero  values 

c.  It  should  have  a  shape  that  varies  from  a 
negative  exponential  to  normal  as  the  mean 
value  of  A',  increases 

d.  The  distribution  of  spacings  between  occur¬ 
rences  should  be  given  by  the  exponential 
disti  bution. 

The  Poisson  distribution  has.  in  fact,  all  these 
characteristics  and  is  given  by 

(Z*  c  11 

/X(V.<»)  -  v,  (4) 

v  =  (>.1.2,3....;  .«>() 
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where  the  parameter  <>  is  the  -.ample  mean,  which 
in  our  ease  varies  I'rom  0.08  lor  lagoons  to  8.0''  tor 
depths  in  excess  of  20  m.  As  we  have  plotted 
,S  10.  these  sample  means  correspond  to  \ 
values  of  0.8  and  80.7  gouges  km  .  Ihe  use  of 
\  10  was  necessitated  by  the  tact  that  \  values 

as  large  as  270  gouges  km  occur.  Ihe  Poisson 
distribution,  on  the  other  hand,  is  not  convenient 
for  values  much  in  excess  of  20.  When  S  10  is 
used,  the  Poisson  probability-  for  an  integer  such 
as  A  is  used  to  represent  the  probability  of  .V  oc¬ 
curring  in  the  interval  25  -s  \  -s  55  gouges  per 
kilometer.  Pxamination  ot  figure  15  shows  that 
the  Poisson  distribution  (the  discrete  values)  does, 
in  tact,  giv  e  a  reasonable  representation  of  the  fre¬ 
quency  plots  of  the  V  values,  although  n  drops 


0  ( nu-st-ei  o*  .]Cuqes  per  lOOml 


oft  too  rapidly  at  large  8  10  values  Ihe  I’ooson 

distribution  also  possesses  the  additive  piopeitv 
that  the  sum  ot  two  independent  Poisson  random 
variables  with  parameters  .<  and  is  also  a 
Poisson  random  variable  with  parameter  - 

<  t  ♦  <  »  ;  . 

Ihe  use  ot  the  Poisson  distribution  brings  to 
mind  us  association  with  the  Poisson  process,  Je 
scribing  the  frequency  ot  i.mdom  events  occurring 
at  a  ciuisiani  rate  al one  a  continuous  space  I  or 
time)  scale.  lo  be  a  Poisson  process,  the  underly¬ 
ing  physical  mechanism  generating  the  events 
must  satistv  the  following  thiee  assumptions: 

1.  Stationaniv  Ihe  probability  ot  at  lea  :  one 
event  in  any  short  interval  is  propo.  oonal  to 
the  length  of  the  interval, 

2.  S'onrmdiipheity — Ihe  probability  of  two  or 
more  events  in  a  short  interval  Av  is  negligi¬ 
ble  in  comparison  to  ui\. 

5.  Independence—  The  number  of  events  in  anv 
interval  is  independent  of  the  number  of 
events  in  any  non-overlapping  interval. 

fhe  probability  distribution  of  the  number  ot 
events  V  in  distance  v  for  a  Poisson  process  is 
given  by 


(51 


n  0.1. 2. 5 .... :  i  v  -0 


w  here  rx  has  replaced  <»  in  cq  4  and  parameter  i  is 
the  average  spatial  rate  of  occurrence  ot  the  even; . 


0  I5r- 


I 


N  ■  i0  (numb®'  cf  gouges  per  ICO  ^  ) 


Figure  If  Relative  frequency  of  different  N  10  values  for  lagoons  and  three  different  water  depth  ranges 
offshore  of  the  harrier  islands  Ihe  solid  line  bar  graph  represents  the  data,  the  discrete  values  indicate  ihe  luted 
Poisson  distribution,  and  the  suppled  area  indicates  the  titled  gamma  distribution. 
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Durance  Be’weer  Gouges  (m; 

f-'mii'V  14.  I  ret/nenev  a!  ucenrrcnce  »'■'  observed  (fciwi'o  between 
the  gouges  oil  Lonely.  IA/vX</. 


W  e  would  judge  that,  when  gouging  is  looked 
on  un  an  annual  event,  it  would  satisfy  the  require¬ 
ments  tor  a  Poisson  process  reasonably  well  as  a 
first  approximation.  We  note  however  that  when 
the  spatial  distribution  of  gouging  is  examined  in 
more  detail  it  is  found  that  there  are  locations 
where  gouges  occur  in  groups  (on  the  seaward 
sides  of  shoals,  for  example).  In  addition,  it  goug¬ 
ing  is  examined  on  a  time  scale  finer  than  yearly, 
the  assumption  of  stationarity  is  clearly  not  satis¬ 
fied,  since  in  many  locations  no  gouging  occurs 
during  the  summer  months.  However,  these  prob¬ 
lems  are  probably  no  worse  than  in  many  other 
areas,  such  as  customer  arrivals  and  number  of 
telephone  calls  per  unit  time,  where  the  Poisson 
process  has  oeen  found  to  be  a  very  useful  model. 

It  is,  of  course,  possible  to  use  other  distribu 
non  functions,  such  as  a  gamma  distribution.  This 
distribution  is  attractive  for  several  reasons,  f  irst, 
it  is  capable  of  assuming  shapes  similar  to  those 
shown  in  f  igure  1.1  (Hahn  and  Shapiro  1967,  Fig¬ 
ures  3 -7b  and  1-8).  It  is  also  an  applicable  distri¬ 
bution  to  data  such  as  V,  that  are  bounded  on  one- 
end  (however  it  is  not  capable  of  treating  the  oc¬ 
currence  of  zero  values).  In  addition,  this  distribu¬ 
tion  has  been  used  successfully  in  a  variety  of  engi¬ 
neering  problems  because  of  its  flexibility  (Benja¬ 
min  and  Cornell  1970).  The  gamma  distribution  is 
given  by 


a  >  o,  \  >  o.  n  >  o 

where  !'(>;)  is  the  gamma  function 

X 

I'd;)  =  j  .v”  1  e  '  (lx.  C) 

u 

Here  the  two  free  parameters  r/  and  X  can  be  con¬ 
sidered  t o  be  shape  and  scale  parameters  respec¬ 
tively.  The  mean,  variance,  and  coefficient  of 
skew  for  the  distribution  are  respectively 


£(.v)  =  >;  X  (8) 

Var(.v)  =  r j/X;  (9) 

1  =  2/ \  ij  .  (10) 

The  exponential  distribution  is,  in  fact,  a  special 
case  of  the  gamma  distribution  with  r  =  1. 


As  can  be  seen  in  Figure  11,  the  gamma  distri¬ 
bution  (the  dashed  lines)  gives  a  very  reasonable 
representation  of  ,V  data  if  the  presence  of  zero 
values  is  arbitrarily  introduced  in  calculating  the 
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I  able  4.  Parameters  of  gamma  distributions  lilted  to  observational  data  on  spatial  gouge  frequency.  lor  the 
combined  offshore  data  set.  the  data  are  expressed  in  terms  of  the  number  ol  gouges  per  kilometer  IP). 


-4: 


appropriate  probabilities.  Note  that  the  gamma 
distribution  is  more  successful  in  luting  the  lareei 
V  10  values  than  is  the  Poisson  distribution, 
which  drops  oil  too  quickly  at  large  values  ol 
\  10.  Table  4  gives  the  values  ol  the  parameters 

ol'  the  titled  gamma  distributions.  I  he  A  and  >/  val¬ 
ues  were  obtained  using  the  maximum  likelihood 
procedure  suggested  by  Thom  (see  Haan  19”.  pp. 


102-106).  In  comparing  the  Poisson  and  the  gam 
ma  mean  values  it  should  be  remembered  that  the 
Poisson  mean  includes  the  etteci'  ol  the  preseme 
ol  zero  \  values  while  the  gamma  mean  does  not. 

I't'mpnral  frequency 

In  investigating  problems  concerning  ice- 
induced  gouging  of  the  sea  floor  it  is  highly  desit- 
able  to  have  independent  information  on  the  rates 


I  able  5.  Number  of  new  gouges  during  the  indieated  time  and  space  intervals.  dtetermined  trnm 

replicate  sonar  data  collected  during  the  summers  of  the  'ears  indicated.* 
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The  symbol  -  indicates  no  data  was  collected. 

Also  given  are  values  ol  a.  the  avetage  number  of  new  gouges  per  kilometer  pet  sear 
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at  which  new  gouges  form  (the  number  ol  new 
gouges  per  kilometer  per  year).  Unfortunately 
such  data  are  rather  limited;  for  our  study  area 
they  are  largely  contained  in  a  paper  h\  Haines  et 
at.  (1978).  Ihis  work  describes  replicate  observa¬ 
tions  made  on  sample  line  35  (see  fig.  5  for  loca¬ 
tion)  during  the  summers  of  1973,  1975,  19'6,  and 
19"  and  on  line  31  during  the  summers  ol  19~5. 
19'6,  and  1977.  We  have  reanalyzed  the  data  set 
from  line  31  for  1976-7'  and  on  line  35  for  the 
1976-77  and  1977_7g  intervals  so  that  the  counts 
of  new  gouges  are  based  on  I -km  sampling  lines. 
We  have  also  analyzed  replicate  runs  on  line  39 
(north  of  Cape  Halkett)  for  1977-78. 

Because  the  quality  of  the  1973  sonar  records 
was  poor  (Reimnitz  et  al.  1977a),  data  based  on 
the  1 97 3 - T’5  time  interval  should  receive  less 
weight  than  the  later  observations.  The  results  of 
this  analysis  and  that  of  Barnes  et  al.  (1978)  are 
combined  and  presented  in  Table  5.  We  have  arbi¬ 
trarily  deleted  the  g  values  obtained  on  line  39  at 
20.3  m  and  farther  offshore,  since  this  portion  of 
the  line  is  known  to  be  in  the  shadow  of  a  nearby 
shoal  area,  thereby  receiving  fewer  gouges.  If  the 
1973-75  data  on  test  line  35  are  also  excluded  be¬ 
cause  of  the  poor  quality  of  the  sonar  record,  we 
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ligure  15.  Values  of  g  (number  of  gouges 
km  vr')  vs  water  depth  (/). 


obtain  an  average  g  value  ol  5.2  gouges  per  kilom 
eter  per  year  with  values  lor  individual  years  varv 
mg  from  2.4  (1975  76)  to  3.5  (1976  77)  to  '  y 
(197'  78).  These  are  appreciably  larger  value' 
than  have  been  obtained  using  similar  procedure' 
oft  the  Mackenzie  Delta  in  15  to  20  m  ol  water 
(0.19  <  0,06  gouges  per  kilometer  per  year.  1  ewis 
19'7a)  and  they  give  a  return  period  per  kilometer 
of  0.2  years  as  compared  to  5.3  years. 

figure  15  shows  a  plot  of  observed  g  values  vs 
water  depth  There  is  no  strong  trend.  In  addition, 
there  is  a  large  scatter  and  zero  values  (I -km  lines 
with  no  new  gouges)  are  rather  evenly  distributed 
at  all  water  depths.  Because  of  this  we  have  treated 
all  the  observations  as  a  single  group. 

f  igure  16  shows  a  plot  of  the  observed  probabil¬ 
ity  of  occurrence  of  different  values  of  g.  The  dis¬ 
tribution  shows  a  strong  positive  skew.  The  Pois¬ 
son  distribution  for  this  set  of  data  is  also  shown. 
The  representation  of  the  data  is  not  encouraging 
(again  the  probability  of  occurrence  falls  oft  much 
too  rapidly  at  large  g  values).  Also  shown  is  a 
gamma  distribution,  which  gives  a  better  fit  (the 
shape  and  scale  parameters  are  respectively  rj  - 
1.205  and  X  =  0.196). 

While  the  characteristics  of  the  new  gouges  are 
being  discussed,  it  is  of  interest  to  examine  the  dis¬ 
tribution  of  their  depths  to  see  it  they  appear  to 
follow  an  exponential  distribution  similar  to  that 
obtained  by  sampling  all  the  gouges  on  the  sea 
floor,  a  data  set  that  contains  a  number  of  old 
gouges  that  presumably  have  been  partially  filled 
with  sediment  as  well  as  new  unfilled  gouges.  The 
observations  used  (n  =  76)  were  from  both  test 
lines  31  and  35  and  occurred  between  1976  and 
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Figure  16.  Relative  frequency  of  different  values 
of  g  (number  of  gouges  km  1  vr  ).  The  discrete  val¬ 
ues  indicate  the  filled  Poisson  distribution  and  the  sup¬ 
pled  area  indicates  the  fitted  gamma  distribution 
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Figure  17.  Semilog  plot  of  relative  fre¬ 
quency  of  occurrence  of  new  gouges  of 
differing  depths  (dj. 


1977.  The  results  are  shown  in  Figure  17.  Again 
the  data  appears  to  show  an  exponential  dropoff 
with  a  a  value  of  4.52  m  '.  This  value  is  dose  to 
but  somewhat  lower  than  the  values  obtained 
from  the  samples  of  all  the  gouges  (taking  1 5  m  as 
a  mean  water  depth  along  the  replicate  sampling 
lines,  we  obtain  a  value  of  5.5  m1  from  Figure  10, 
as  contrasted  with  4.5  nr'  from  the  new  gouges). 
That  new  gouges  should  have  a  lower  X  value  than 
a  corresponding  distribution  of  old  and  new 
gouges  could  be  anticipated  (E.  Phifer,  pers. 
comm.)  from  the  observation  that  at  other  loca¬ 
tions  deep  gouges  in  the  sea  floor  receive  more  fill 
per  year  than  do  shallow  gouges  (Fredsoe  1979). 
At  the  present  there  clearly  is  no  strong  reason  to 
doubt  that  the  distribution  of  new  gouge  depths  is 
exponential  or  that  the  X  values  that  will  be  ob¬ 
tained  are  greatly  different  (presumably  slightly 
less)  than  values  obtained  from  our  earlier  analysis 
of  all  the  gouges. 


Extreme  value  analysis 

Another  way  to  view  portions  T  the  gouging 
data  is  by  extreme  value  analysis.  In  this  case  the 
complete  data  set  is  not  examined.  Instead  the 
largest  (or  smallest)  value  in  each  of  a  number  of 
specified  sampling  intervals  is  used.  In  most  appli¬ 
cations,  such  as  hydrology,  the  data  are  in  the 
form  jf  time  series,  and  the  largest  (smallest) 
event  in  each  of  a  sequence  of  fixed  time  inters  als 
is  used  to  generate  a  distribution  of  rare  events.  In 
our  study,  the  basic  data  set  is  a  space  series,  as 
separate  frequency  distributions  of  gouge  charac¬ 
teristics  were  developed  for  each  kilometer  of 
sampling  line.  For  instance,  in  a  kilometer  of  line 
one  might  observe  85  gouges  of  different  depths, 
with  the  largest  gouge  having  a  value  of  2.2  m;  in 
the  next  kilometer  there  might  be  178  gouges  with 
a  maximum  value  of  3.1  m.  The  extreme  value 
distribution  would  then  be  composed  of  the  values 
2.2,  3.1,  and  subsequent  values.  Good  discussions 
of  the  different  types  of  extreme  value  distribu¬ 
tions  can  be  found  in  Hahn  and  Shapiro  (1967), 
Benjamin  and  Cornell  (1970)  aid  Haan  (1977). 

The  particular  extreme  valui  distribution  ap¬ 
plicable  to  a  given  situation  depends  on  the  nature 
of  the  initial  distribution  being  sampled  and  on  the 
sample  si/e  n,  with  the  extreme  distribution  being 
approached  asymptotically  as  n  becomes  large.  A 
common  problem  is  that  many  times  n  does  not 
appear  to  have  been  large  enough,  and  the  extreme 
value  distribution  that  would  be  expected  to  apply 
to  a  given  data  set  is  not  particularly  successful  in 
fitting  it.  For  instance,  a  Type  I  extreme  value  dis¬ 
tribution  should  apply  to  maximum  values  sam¬ 
pled  from  an  initial  distribution  that  is  of  the  ex¬ 
ponential  type.  However,  Tucker  et  al.  (1979),  in 
their  study  of  maximum  pressure  ridge  heights 
whose  initial  distribution  appears  to  be  the  expo¬ 
nential  type,  found  that  their  data  were  not  linear 
on  Type  I  paper  but  were  effectively  linearized  by 
standard  probability  paper.  Similar  results  have 
been  obtained  by  other  workers  in  hydrology  and 
in  Monte  Carlo  simulations  by  Slack  et  al.  (1975). 
In  practice,  a  number  of  different  approaches 
(Type  I,  normal,  log-normal,  log-Pearson  Type 
Ml)  are  commonly  tested,  and  the  most  successful 
relation  is  selected  to  analyze  the  data. 

Gouge  depths 

As  we  have  shown,  gouge  depths  appear  to  be 
exponentially  distributed.  Therefore,  the  appro¬ 
priate  extreme  value  distribution  for  maximum 
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Figure  IS.  Plots  of  dmav  v.v  water  depth  (?)  for  different  regions  within  the  study  area. 


gouge  depths  should  be  a  Type  I  distribution. 
However,  testing  shows  that  the  data  were  not  lin¬ 
earized  by  either  a  Type  I,  a  normal,  or  a  log-nor¬ 
mal  distribution.  However,  a  log-Pearson  Type  111 
(LPIII)  distribution  proved  to  be  quite  effective. 
This  distribution,  which  is  in  fact  a  three- 
parameter  gamma  distribution  fitted  to  the  log,0 
of  the  extreme  values,  has  been  used  successfully 
in  treating  flood  observations  (USWRC  1977). 
The  three  parameters  describing  an  LPIII  distribu¬ 
tion  are  the  mean  A,  the  standard  deviation  S,  and 
the  skew  coefficient  G,  which,  if  X  s  log,,  rfmax, 
where  dmax  is  the  maximum  gouge  depth  in  a  kil¬ 
ometer  track  and  N  is  the  number  of  maximum 
gouges,  are  calculated  as  follows: 


c  UX-X)! 

*  ~  N- 1 

NL(X-X)' 

(N-  l)(N-2)S' 

The  computed  dmax  value  is  then  given  by  the  rela¬ 
tion 

\og«,dmax  =  X-t  KS  (II) 

where  K  is  the  Pearson  Type  III  coordinate  ex¬ 
pressed  in  magnitudes  of  the  standard  deviation 
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from  the  mean  for  various  exceedance  percent¬ 
ages.  Values  of  K  are  functions  of  C  and  are  given 
in  Appendix  3  in  USWRC  (1977),  as  are  the  com 
puling  equations  for  \  .S',  and  O'. 

In  analyzing  the  r/max  values  on  gouging,  indi¬ 
vidual  plots  (Fig.  18)  were  prepared  show  ing  </1)1JX 
vs  ;  for  five  different  areas.  The  different  regions 
were  compared  by  overlaying  the  figures  on  a  light 
table.  If  differences  in  water  depth  are  taken  into 
consideration,  the  data  from  l  onely.  Harrison 
Bay.  Jones  Islands,  and  McClure  Islands  overlap 
very  well  and  appear  to  form  one  continuous  dis¬ 
tribution.  Therefore,  as  before,  the  data  were 
pooled  into  one  sample.  The  data  from  the 
lagoons  were  treated  separately,  both  because  they 
appear  different  and  they  represent  a  different 
marine  environment. 

Another  characteristic  of  the  </max  data,  which 
might  be  anticipated  from  our  earlier  discussion 
and  that  is  apparent  in  Figure  18,  is  that  the  values 
clearly  change  with  water  depth.  There  are  null 
values,  many  small  values,  and  no  large  values  in 
shallow  water;  large  values  of  r/max  become  in¬ 
creasingly  common  with  increasing  ;;  and  small 
values  are  rare  in  water  deeper  than  20  m. 
Therefore,  as  before  the  pooled  offshore  </max 
data  were  separated  into  5-m  water  depth  in¬ 
crements.  As  no  similar  trend  was  apparent  in 
the  data  from  lagoons  and,  as  the  depth  range  was 
limited,  these  results  were  not  separated  into 
similar  groups. 

In  analyzing  the  data,  two  problems  were  en¬ 
countered.  First,  in  a  number  of  shallow  water 
areas  we  commonly  found  appreciable  lengths  of 
track  that  did  not  contain  gauges,  resulting  in 
(fmax  =  0  values.  For  instance,  in  the  data  set  for 
lagoons,  119  km  of  the  324  km  sampled  (37#’<i) 
were  gouge-free.  This  precludes  the  normal  statis¬ 
tical  analysis  of  the  data  using  an  LPIII  distribu¬ 
tion,  as  the  log, i,  of  zero  is  minus  infinity.  Sec¬ 
ondly,  in  a  number  of  cases  it  was  impossible  to 


determine  precisely  the  depth  o!  the  miuIici 
gouges,  only  that  a  gouge  existed  and  that  it-, 
depth  was  less  than  some  specified  value  Stub 
gouge  depths  are  identified  by  circles  in  I  iguie  18 
In  most  cases  they  had  values  ot  less  than  IU  in 
and  were  situated  in  shallow  water.  This  created 
considerable  uncertainty  in  specifying  the  exact 
number  ot  gouges  in  the  (I  I-  and  0  2-m  depth 
classes.  Where  such  gouges  were  common  tat 
water  depths  of  less  than  10  ml,  large  ( »  values  and 
1  I’ll!  distributions  were  obtained  that  were  not 
particularly  good  fits  to  the  data  at  the  larger  r/n 
values  (which,  of  com  se.  ,s  the  area  of  prime  inter¬ 
est). 

Both  of  these  problems  were  handled  using  a 
procedure  developed  for  treating  zero  flood  years 
and  incomplete  records  in  hydrology.  First,  the  (I. 
0.1-,  and  l).2-m  values  were  deleted  from  the  sam¬ 
ples.  Then  the  A.  S',  and  (/  parameters  were  calcu¬ 
lated  from  the  censored  distributions  and  used  to 
calculate  </  as  a  function  of  exceedance  proba¬ 
bility.  These  exceedance  probabilities  were  then 
adjusted  by  multiplying  them  by  the  ratio  ot  the 
number  of  values  in  the  censored  distribution  to 
the  number  of  values  in  the  uncensored  dixtribu 
lion  (i.c.  with  the  0.  0.1.  and  0.2  values  included). 
The  results  were  then  plotted  on  log-probability 
paper  for  comparisons  with  the  observed  data.  In 
plotting  the  data  against  the  adjusted  curve,  the 
plotting  positions  were  determined  bv  using  the 
Weibull  plotting  formula 

P  =  m  t.V+  1) 

where  P  is  the  exceedance  probability;  m  the  se¬ 
quence  of  rf  values,  with  the  largest  values  cor¬ 
responding  to  in  =  1,  the  next  largest  value  corre¬ 
sponding  to  in  =  2,  etc.;  and  ,V  the  total  number 
of  data  points  before  censoring  (i.c.  including  0. 
0. 1 ,  and  0.2  values). 


Table  6.  Parameters  of  the  log-Pearson  Type  III  distribution  (determined  from  values  of  </milx  observed  along 
l-Wm  sampling  lines.  Values  outside  (he  harrier  islands  inetude  data  trom  Harrison  Bav  and  north  of  l  onelv. 
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Table  6  gives  t he  A.  .S',  and  (/  values  calculated 
from  the  different  sets  of  censored  data  as  well  as 
the  adjustment  ratio  and  the  number  of  < im  val¬ 
ues  equal  to  zero  and  between  0.3  and  zero.  The 
exceedance  probabilities— the  probabilities  that 
given  a  single  kilometer  of  sample  track,  the  maxi¬ 
mum  gouge  depth  will  be  equal  to  or  greater  than 
some  indicated  value,  r/max— are  shown  in  figure 
19.  Also  shown  is  the  spatial  recurrence  interval 
for  1-km  segments  with  one  or  more  exceedances, 
which  is  equal  to  the  reciprocal  of  the  exceedance 
probability.  This  parameter  gives  the  expected 
number  of  kilometers  of  sea  floor  that  must  be  ob¬ 
served  before  the  maximum  gouge  depth  in  one  of 
(hose  kilometers  is  expected  to  equal  or  exceed 
z/max.  Another  parameter  of  possible  interest  is  the 
number  of  kilometers  per  100  km  of  sample  track 
in  which  the  maximum  gouge  depth  is  expected  to 
equal  or  exceed  r/max.  This  number  can  be  ob¬ 
tained  simply  by  multiplying  the  appropriate  ex¬ 
ceedance  probability  by  100.  The  curves  sweep 
across  the  graph  and  show  systematic  changes 
with  water  depth  as  was  expected.  The  10-  to  15- 
and  the  20-  to  25-m  curves,  which  arc  omitted  to 


restrict  clutter,  lie  as  expected  on  the  figure.  I  he 
30-  to  35-m  curve  is  very  similar  to  the  25-  to  30-m 
curve,  which  is  not  too  surprising  as  there  are  not 
mans  </mjx  values  in  the  30-  to  35-m  range.  It 
should  also  be  noted  that  in  the  plots  of  </(llax  vs 
from  the  Mackenzie  Delta  region  (I  ewjs  )9"a), 
the  </m  values  peak  at  approximately  40  m  and 
decrease  in  deeper  water. 

In  f  igure  19  the  0  to  5-m  data  and  the  data  from 
the  lagoons  overlap  each  other.  As  there  are  only 
three  data  points  in  the  0  to  5-m  data  set  (as  the 
result  of  censoring  the  lower  values),  the  calcu¬ 
lated  curve  was  not  particularly  similar  to  the 
curves  front  deeper  water.  The  curve  presented  in 
f  igure  19  is  based  on  the  data  from  lagoons  and 
appears  to  give  a  reasonable  representation  ol  the 
0  to  5-m  data  points  as  well. 

f  igure  20  presents  .V  /o.ez/m.lx,  -I.  O.  and  .5 
plotted  as  a  function  of  c.  This  plot  should  be  use¬ 
ful  to  those  interested  in  developing  eq  1 1  to  applv 
to  other  water  depth  intervals  than  those  con¬ 
sidered  here.  The  most  systematic  change  in  a  pa 
rameter  with  c  is  the  roughly  linear  increase  in  A 
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Inure  20.  Parameters  relative  to  the  determination 
of  eq  II  shown  as  a  function  of  water  depth  (//. 


JC 


figure  21  Plot  of  w  (  for  I -km  line  segments  r% 
water  depth  0 j  for  all  locutions  except  lagoons. 
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Gouge  widths 

Figure  21  shows  all  maximum  gouge  widths 
< w max)  nieasurt'd  outside  the  barrier  islands  com¬ 
pared  with  the  average  water  depth  (c).  The  trends 
are  similar  to  those  present  in  Figure  18,  which 
plotted  dnuiy  vs  average  water  depth.  There  is  a 
general  increase  in  «m.iv  as  ;  increases.  This  may 
simply  reflect  that,  on  the  whole,  gouges  that  are 
deeper  are  also  wider.  In  addition,  in  deeper  water 
there  do  not  appear  to  be  any  small  u-  x  values  as 
there  were  in  shallow  water. 

Lateral  embankment  heights 

Finally,  a  comparison  of  /i  ,  the  maximum 
lateral  embankment  heigftt,  and  </m  is  presented 
in  f  igure  22  (the  numbers  indicate  the  number  of 
values).  It  is  hardly  surprising  that,  on  the  aver¬ 
age,  regions  with  deeper  gouges  should  contain 
higher  embankments  as  the  material  front  the 
gouges  produces  the  embankments.  However,  we 
were  surprised  at  how  symmetrically  the  values 
were  distributed  around  the  one-to-one  line.  This 
is  show  n  by  the  histogram  (see  the  inset  in  Fig.  22) 
of  the  relative  frequency  of  deviations  from  the 
one-to-one  line  (measured  normal  to  that  line). 


APPLICATIONS  TO 
OFFSHORE  DESIGN 

In  the  preceeding  sections  we  have  attempted  to 
systematize  and  clarify  some  of  the  essential  char¬ 
acteristics  of  a  large  set  of  measurements  on  the 
geometry  of  ice-induced  gouges  in  the  sediments 
of  the  Alaskan  portion  of  the  shelf  of  the  Beaufort 
Sea.  These  observations  are,  of  course,  valuable  in 
themselves.  For  instance,  it  is  useful  to  know  that 
outside  of  the  barrier  islands  in  water  up  to  38  m 
deep  the  deepest  gouge  observed  was  2.6  tn,  ob¬ 
tained  from  a  sample  of  20,313  gouges  collected 
over  1500  km  of  sampling  track.  In  the  protected 
lagoons,  on  the  other  hand,  the  deepest  gouge  (0.7 
m)  was  much  shallower  (from  a  sample  of  41 
gouges  obtained  from  298  km  of  sainplint  'rack) 
and  a  large  percentage  (92®7o)  of  the  l -km  ,eg 
ments  examined  contained  no  gouges  at  an.  In  the 
remainder  of  this  section  we  will  attempt  to  use  the 
data  analysis  performed  earlier  in  this  paper  to 
make  a  scries  of  preliminary  estimates  of  the  prob¬ 
ability  of  occurrence  of  gouges  with  certain  pre¬ 
scribed  depths  and  frequencies. 

Gouge  depths 

To  obtain  the  exceedance  probability  for  the  oc¬ 
currence  of  gouges  of  different  depths,  given  that 
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L inure  J.L  Plot  oj  the  exceedance  probability  /t >  \t\i/ 
vs  volute  depth  for  different  water  depths  (/ 1  in  the 
offshore  region  unprotected  In  hat  tier  islands,  b:  n// 
cidannit  the  Uvfxy  values,  the  \  values  arc  obtained 
from  the  relation  shown  in  liiture  1(1. 


gouging  has  occurred,  the  relation  in  Figure  JO  can 
be  used  to  obtain  an  estimate  of  X  applicable  to  the 
water  depth  of  interest.  The  exceedance  probabil¬ 
ity  is  then  obtained  from  eq  2.  For  instance,  lor  a 

A 

water  depth  of  5  m,  X  -  8.16  and 

P\D>  I)  =  exp[  8.16(1-0.2))  146-10' 

gives  the  probability  of  a  gouge  exceeding  I  m  in 
depth.  Therefore,  using  eq  3,  one  gouge  in  685 
would  be  expected  to  be  at  least  I  in  deep.  The 
0. 2-m  correction  in  the  above  calculation  is  caused 
by  the  fact  that  (he  0  to  0  2-m  depth  class  was  dc 
leted  in  the  estimation  of  X.  At  the  same  water 
depth,  one  gouge  in  2.39  million  would  be  expect¬ 
ed  to  be  at  least  2  m  deep.  For  35  m  of  water  (X 
2.46)  things  are  very  dif  ferent;  one  gouge  in  seven 
exceeds  1  m  and  one  in  980  exceeds  3  m.  A  graphic 
display  of  the  variations  in  the  exceedance  proba 
bility  as  a  function  of  water  depth  for  the  offshore 
region  is  given  in  Figure  23. 


The  \  values  determined  tor  lagoons  appear  to  gouge  greater  than  or  equal  to  i /  Nevertheless, 

he  in  the  7-  to  9-m  range,  i.e.  in  general  agree-  both  approaches  itin  be  useful  it  applied  apptopi  i- 

nient  with  the  X  values  obtained  from  similar  ately. 

water  depths  in  the  offshore  data  set.  C  onsider  three  20  km  pipeline  routes,  one  in  the 

lagoons  and  two  at  sites  unprotected  In  islands  in 
F.xlreme  value  statisties  5  to  It)  and  I*)  to  20  m  ot  water  respectively.  1  or 

It  is  important  to  note  two  factors  concerning  the  lagoons,  the  extreme  value  exceedance  proba¬ 
tive  extreme  value  statistics  that  have  been  pre-  brlitv  tor  1  km  sampling  intervals  is  approximately 

seated,  first,  the  sampling  lines  cross  the  gouges  0.0065  and  0  00012  toi  gouge  depths  ot  0.5  and 

at  a  variety  of  angles.  Therefore,  from  an  area  1 .0  m  respectively .  corresponding  to  spatial  recur - 

where  the  gouging  is  strongly  aligned,  the  maxi-  mice  intervals  ot  154  and  ’6') 2  km.  C  orrespond 

mum  value  used  was  selected,  in  some  cases,  from  ing  values  for  5  to  10  in  and  2'  ro  lO-m  water 

a  small  number  of  gouges  (when  the  sampling  line  depths  outside  ot  the  barrier  islands  are  given  in 

nearly  paralleled  the  gouges)  and  in  other  cases  Table  ".  Based  on  this  table  we  could  conclude 

from  a  much  larger  number  (when  the  sampling  that  if  one  was  to  contemplate  using  an  engineer- 

was  perpendicular  to  the  gouges).  We  have  not  at-  ing  technique  that  would  encounter  difficulties  m 

tempted  to  correct  the  extreme  value  data  in  the  the  presence  of  gouges  ol  1  m  or  more,  we  would 

manner  that  we  corrected  the  observations  on  the  not  anticipate  problems  in  constructing  a  20-km 

observed  number  of  gouges  per  kilometer  (  Y)  to  line  within  the  lagoons.  On  the  othci  hand,  at 

the  number  that  would  be  expected  if  the  sampling  water  depths  ol  25  to  2()  m  we  would  expect  to  en 

were  perpendicular  to  the  gouging  (  V  ).  We  do  not  counter  gouges  at  least  I  m  deep  in  roughly  15  ot 

know  how  to  make  such  a  correction.  the  20  km. 

Secondly,  it  should  be  realized  that  the  extreme  Another  parameter  ol  inteiest  is  the  probability 

value  and  the  complete  distribution  techniques  /'(.  I)  that  the  maximum  gouge  depth  per  kilometei 

give  estimates  of  two  different  factors.  Die  ex  will  equal  or  exceed  a  given  value  (e.g.  I  ml  along 

treme  value  approach  provides  an  estimate  of  the  the  pipeline.  This  is  calculated  as  follows: 

number  of  I -km  segments  that  will  have  at  least 

one  gouge  greater  than  or  equal  to  some  specified  />(  l>  1  /*(/() 

value  r/maK  along  a  given  length  of  sampling  line 

On  the  other  hand,  an  estimate  using  the  complete  where  !'(H)  is  the  prohabihtv  t hat  the  maximum 

PDF  gives  the  expected  number  of  gouges  along  gouge  depth  per  kilometer  will  not  equal  or  exceed 

the  line  that  are  greater  than  or  equal  to  r/  .  The  I  m  in  any  of  the  20  km  /’(H)  in  tuin  equal'  the 

two  estimates  arc  not  the  same  because  a  given  probability  that  the  maximum  gouge  depth  pet  ki 

1-km  sampling  segment  may  have  more  than  one  lometer  will  not  he  cl  m  in  the  lust  kilomcict 


Table  7.  I.xceedance  probabilities.* 
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*(ti\cn  I  km  ol  sample  track.  spatial  recurrence  intervals  tor  )  km  segments,  and 
probabilities  /VOrhar  the  maximum  gouge  depth  per  kilometer  will  equal  ot  ex 
coed  the  mdieitlcd  gouge  depth  along  a  20  km  line  based  on  the  extreme  value 
statistics 
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multiplied  in  the  ptohubilitv  that  it  mil  no!  he  •  ! 
m  itt  (he  second  kilometer,  ete.  up  to  the  twentieth 
kilonietet.  Assuming  that  eaeh  kilomctci  has  the 
same  prohahilit>  /’(O  that  the  ni.iMiniiiii  gouge 
depth  per  kilometer  will  not  he  -1  m,  then  l\H) 
I/’IO)  )\l  )  is,  however.  equal  to  I  minus  the 
probability  /’(/))  that  the  maximum  troupe  depth 
pet  kilometer  will  he  '1  m.  In  short. 

l\  1)  I  (I  /’(/>)!"  (12) 

where  n  is  the  number  of  I -km  segments  compos- 
mg  the  line.  In  our  example,  n  20  and  l\P) 
0.0001 3  for  lagoons  so  that  /’(.  l)  0. 002(1.  These 

values  as  well  as  similar  values  at  watet  depths  of  5 
to  10  m  and  25  to  50  in  are  also  included  in  fable 

As  is  shown,  the  probability  of  encountering  an 
extreme  gouge  with  a  depth  equal  to  or  greater 
than  4  m  in  water  25  to  30  in  deep  is  appreciably 
larger  than  the  probability  of  encountering  a  1-m 
extreme  gouge  in  the  lagoons. 

Burial  depths 

fhe  problem  of  burial  depth  can  be  considered 
m  several  dilferent  ways.  Clearly,  every  gouge 
greater  than  a  specified  value  is  important,  so  it  is 
necessary  to  use  the  I’D!  based  on  the  complete 
set  of  gouge  depths,  as  opposed  to  the  extreme  v  al 
ue  distribution  based  on  the  maximum  gouge  in 
each  kilometer. 

f  irst,  we  consider  the  situation  where  we  wish  to 
bury  the  pipeline  at  a  depth  so  that  it  is  all  covered 
(assuming  an  acceptably  low  probability  of  en¬ 
countering  a  gouge  deeper  than  out  burial  depth 
that  would  leave  the  line  uncovered).  In  this  ease 
we  are  dealing  with  gouge  depths  as  they  exist  on 
the  sea  floor  at  a  given  instance  of  lime.  Again  as 
an  example  we  will  consider  a  20  km  line  that  will 
be.  in  turn,  restricted  to  lagoons  and  to  water 


depths  of  5  to  10  m  and  25  to  30  in  outside  the  bat 
iiei  islands.  We  will  also  consider  the  case  wheie 
the  diiection  ol  the  line  is  20  off  the  direction  of 
the  gouges  as  well  as  nottnal  to  the  diiection  ol  the 
gouges,  foi  instance,  al  a  water  depth  of  2^  to  30 
ill  we  would  expect  to  encounter  an  average  of  SO 
gouges  per  kilometer  if  the  line  is  normal  to  the 
gouges  and  SO  siti  20  27  gouges  pci  kilonietet  0 

the  angle  between  the  gouges  and  the  hue  is  20  . 
Considering  20- km  lines,  this  corresponds  to  lhOO 
and  540  gouges,  tespeclivelv  Next,  vine  must  de¬ 
cide  how  many  gouges  can  be  inlet ated  deepen 
than  the  depth  of  burial.  We  will  take  two  cases: 
one  exceedance  pet  20  km  and  one  exceedance  per 
100  km.  Burial  depths  (\)  can  then  be  calculated 
from  eq  3.  which,  when  real  ranged  and  modified 
tv'  neat  the  above  cases,  becomes 

u(/).s,/i  h|/)>«/J  . 

A  V  (sin  H)l  v 

or,  rearranging 

v  0.2  '  In 


ill/)  -</] 
V  (sin  f/)/ 


(14) 


As  stated,  at  a  water  depth  of  5  to  10  m,  X  '.3, 
V  10,  0  20  or  90  .  I  20  or  100  km  and 

/;|/)'</}  I  inasmuch  as  we  only  wish  to  allow 
one  exceedance.  The  results  of  several  such 
calculations  are  given  in  fable  8. 

Unfortunately,  the  problem  we  would  really  like 
tv'  solve  is  somewhat  different  and  more  difficult 
than  the  above;  a  pipeline  is  buried  and  we  wish  to 
estimate  as  a  function  of  burial  depth  how  often 
(in  a  time  sense)  the  pipeline  can  be  expected  tv'  be 
impacted  by  a  pressure  ridge  keel.  Ibis  problem 
also  requires  knowledge  of  the  rates  of  occurrence 


Table  8.  Kstimated  burial  depths  assuming  that  one  existing  gouge  will  exceed  the  burial 
depth  along  the  length  of  the  line. 
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of  new  gouges.  What  length  of  time  chi  the  ob¬ 
served  gouge  sets  represent'.’  i  his  question  van  be 
examined  from  several  different  viewpoints.  I  nst, 
we  can  estimate  sedimentation  tales  m  the  studs 
area  to  see  how  last  gouges  would  he  erased 
(filled),  assuming  uniform  sedimentation.  Avei 
age  sedimentation  rates  appear  to  he  unite  low 
Keimnu/  et  at.  ( I*)"”)  obtained  an  average  value  ot 
0.06  cm  vr  by  dividing  the  observed  average 
thickness  of  recent  (Holocene)  sediments  (3  m>  hv 
the  period  ot  time  then  study  area  was  believed  to 
have  been  covered  hv  the  sea  (5000  years).  1  evns 
(1977a)  obtained  similar  hut  generally  highei  val 
ties  (0.05  to  0.2  cm  yr  )  for  his  study  area  north  of 
the  Mackenzie  Delta.  I 'sing  the  006  cm  yt  value 
and  assuming  that  no  other  processes  are  active,  it 
would  lake  about  1666  years  to  fill  a  I -in  deep 
gouge  and  5000  years  to  fill  a  3- m  deep  gouge 
Based  only  on  this  information,  an  observed 
gouge  set  would  represent  a  long  petiod  ot  time 

In  the  above,  the  assumption  ot  unifotm  sedi¬ 
mentation  v'n  the  shell  is  probably  in  emu  A 
gouged  bottom  morphology  cteates  abrupt  local 
reliet  and  local  sedimentation  t ate  anomalies  that 
amount  tv'  large  differences  in  sedimentation  rates 
over  short  distances.  Gouge  embankments  may  be 
sites  of  erosion  while  the  gouges,  as  depressions, 
act  as  loci  of  much  higher  rates  of  sedimentation 
than  would  be  apparent  on  a  regional  basis,  fur¬ 
thermore,  sedimentary  structures  tn-shore  of  20  m 
show  shelf  deposits  to  consist  of  gouge  infill  mate¬ 
rial  (Barnes  and  Reimnitz  1*174.  Barnes  et  al. 
1979). 

In  addition,  it  is  becoming  increasingly  ap¬ 
parent  that  shallow  water  gouges  are  rapidly  oblit¬ 
erated  due  to  high  levels  of  hydrodynamic  activity 
(Kovacs  1*172.  I’ilkington  and  Marcellus  19X1). 
f  or  instance,  recent  field  observations  by  Barnes 
and  Reimnitz  (1979)  show  that  the  extensive  open- 
water  conditions  that  occurred  during  the  summer 
of  1977  resulted  in  hydrodynamic  conditions  (pre¬ 
sumably.  large  waves  and  wind-generated  shelf 
currents  associated  with  the  presence  of  a  large 
fetch)  that  have  obliterated  ice  gouges  to  a  water 
depth  of  13  m  and  caused  pronounced  infilling  of 
gouges  in  deeper  water.  Apparently,  the  rates  of 
reworking  a. id  redepositing  sediment  from  such 
episodic  events  are  much  greater  than  the  average 
sediment  accumulation  rate  on  the  Beaufort  Sea 
shelf. 

We  know  of  no  studies  of  the  recurrence  fre¬ 
quency  of  conditions  such  ar  those  observed  dur¬ 
ing  the  summer  of  1977.  but  we  would  guess  that 
they  are  fairly  common,  with  return  periods  of  no 


moie  than  25  yeats.  I  wenty-live  years  appears  to 
be  a  reasonable  estimate  for  the  return  period  ot 
significant  stotm  surges  along  the  coast  of  the 
Beaufort  Sea  (Reimnitz  and  Maurer  197X),  events 
that  would  presumably  be  associated  with  similar 
or  more  energetic  hydrodynamic  conditions  In 
short,  although  sedimentation  rates  might  lead 
one  to  believe  that  the  Beaufort  Shelf  is  a  rather 
static  environment  sedimentologically .  this  is  lar 
from  the  case;  particularly  in  locations  where 
water  depths  are  less  than  10  m. 

I  hetelore.  m  most  of  the  area  we  have  studied, 
we  would  not  have  confidence  in  the  assumption 
that  the  sea  floor,  as  seen  at  a  given  time,  repre¬ 
sents  a  steady  state  condition  with  the  number  of 
new  gouges  per  unit  time  equaling  the  number  ot 
gouges  infilled  by  sedimentation  plus  the  number 
of  new  gouges  superimposed  on  existing  gouges. 
Such  statistical  time  invariance  ol  the  gouging  i' 
an  essential  assumption  it  the  rate  of  production 
v't  new  gouges  is  estimated  using  the  scour  budget 
approach  developed  by  1  ewis  (197“a.h).  \\  e  think 
the  method  is  interesting  and  quite  possibly  appli¬ 
cable  to  certain  regions  v't  gouging,  for  instance 
ottshore  areas  in  the  Chukchi  Sea  with  water 
depths  ot  30  tvi  50  m.  However,  tot  the  Bvautori 
Sea  in  general,  and  in  particular  for  water  depths 
less  than  20  m.  hi  fee)  that  the  applicability  ol  the 
steady-state  assumption  is  doubtful. 

Another  approach  used  to  get  a  rough  estimate 
of  the  age  of  an  observed  set  of  gouges  is  to  divide 
the  average  value  for  the  annual  sum  of  the  gouge 
widths  by  the  length  of  the  sample  track  (Reimnitz 
et  al.  197'a).  For  instance,  if  our  sample  line  is  10 
km  long  and  we  obtain  an  average  of  500  m  ot  new 
gouges  crossing  the  line  each  year,  we  then  take  20 
years  as  an  estimate  of  the  time  period  in  which 
the  gouges  are  completely  replaced.  In  fact,  such 
estimates  give  the  shortest  period  of  time  in  which 
the  gouge  set  could  be  replaced  (an  event  of  very 
low  probability),  as  ice  presumably  plows  the  sea 
floor  in  a  random  (not  a  systematic)  manner. 
Therefore,  the  fact  that  a  given  segment  of  a  line 
has  just  been  gouged  has  no  effect  on  the  proba¬ 
bility  that  the  segment  will  be  gouged  the  next  year 
(or  the  next  month). 

Still  another  approach  using  the  same  data  set 
assumes  that  an  increasingly  large  proportion  of 
the  bottom  is  regouged  before  the  entire  bottom  is 
gouged  (Barnes  et  al.  1978).  In  this  scheme,  if  10°« 
of  the  seabed  is  gouged  each  year  then  in  the  first 
year  10°o  is  impacted  with  new  gouges  but  in  the 
second  year  only  1 9(ro  is  gouged  as  1  °o  of  the 
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gouges  occurred  in  areas  already  gouged.  I  his  can 
be  expressed  as  the  polynomial 

(>,  1  (l  A)7  (15) 

where  Ot  is  the  traction  of  the  bottom  gouged 
since  7„,  A  is  the  fraction  of  the  bottom  gouged 
each  year,  and  /'is  the  time  in  years  measured  rel¬ 
ative  to  T„. 

Finally,  attempts  have  been  made  to  combine 
information  on  pressure  ridge  keels,  pack  ice  drift, 
and  observed  distributions  of  gouge  depths  to  esti¬ 
mate  required  burial  depths  (Pilkington  and  Mar- 
cellus  1981.  Wadhams  1983).  As  the  first  two  of 
these  parameters  are  very  poorly  known,  such  esti¬ 
mates  are  uncertain.  This  technique  also  appears 
to  give  maximum  gouge  depths  that  are  apprecia¬ 
bly  deeper  than  observed.  More  will  be  said  about 
this  later. 

We  believe  that  at  present  to  examine  adequate¬ 
ly  the  problem  of  pipeline  burial,  independent  in¬ 
formation  on  gouging  rates  and  the  depths  of  re¬ 
cent  gouges  is  essential.  As  we  have  described,  our 
information  on  this  subject  is  hardly  what  we 
would  desire.  Nevertheless  it  is  enough  to  allow  us 
to  make  an  initial  approach  to  estimating  burial 
depths.  To  summarize  our  observations  on  recent 
gouges,  we  found  that  g,  the  number  of  gouges  per 
km  per  year,  varied  from  2.4  to  7.9  with  a  mean  of 
5.2.  There  also  was  no  apparent  relation  between  g 
and  water  depth.  The  PDF’  for  recent  gouges  was 
exponential  with  a  X  value  of  4.5  m  ',  a  value  that 
is  I  m  1  less  than  comparable  X  values  from  all  the 
gouges  existing  on  the  sea  floor  at  a  given  time. 


I  sing  this  into!  nation  we  can  now  make  pic 
limmaty  estimates  ol  t tie  biin.il  depths  reqimcd  m> 
that  a  pipeline  ol  a  eiven  length  will,  on  llic  avei 
age.  be  gouged  once  dining  a  specilied  penod  ol 
time  1 1 oi  iiistaiuc.  one  time  in  IOO  or  in  It  too 
years).  1  o  do  this.  Inst  estimate  X .  the  total  mini 
bet  ol  gouges  that  will  occui  dining  the  pioposcd 
lifetime  eit  the  pipeline  b\ 


\  g  /  I  sin" 


(16) 


where  g 

/ 

1 

0 


average  number  ot  gouges  km  >r 
occurring  along  the  pipeline  route 
proposed  lifetime  in  seals 
length  of  the  line  in  kilometeis 
angle  bet  wen  the  route  and  the  trend 
of  the  gouges. 


As  we  only  consider  one  contact  in  /.  V|/)><7] 
in  eq  3  equals  1  and  we  obtain 
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g  T l  sin 6 


or 
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(18) 


In  Table  9  we  show  a  series  of  burial  depth  esti¬ 
mates  made  using  eq  18.  In  these  calculations  we 
have  used  both  the  observed  *X  value  for  the  exist- 
ing  gouge  set  from  Figure  10  as  well  as  X  1  as  an 
estimate  of  the  corresponding  parameter  for  new 


Table  9.  Kstimaled  burial  depths  assuming  one  contact  between  a  pressure  ridge  keel  and  the  pipeline  in  100 
years.  It  ulculltions  made  using  eq  18.1 
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Table  HI.  (  ompurisons  between  burial  iteplhs  In  (he  lop  of  a  76-km  pipeline  for  a 
1000-,  100-,  and  10-year  relurn  period.* 
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gouges.  In  using  the  table,  note  that  a  20-year  life¬ 
time  for  a  100-km  line  is  identical  to  a  100-year 
lifetime  for  a  20-km  lit"-  As  can  be  seen  in  the 
table,  it  is  very  import,!  10  obtain  data  that  will 
allow  improved  estimates  of  X  and  g  for  new 
gouges.  In  general  it  can  be  said  that  slight  in¬ 
creases  in  the  burial  depth  (a  few  tens  of  centime¬ 
ters)  result  in  appreeiable  increases  in  the  safety  of 
the  line.  This  statement  is  particularly  true  in  shal¬ 
low  water  where  X  is  large. 

In  Table  10  we  have  also  included  a  comparison 
between  our  estimates  of  burial  depths  and  those 
of  Pilkington  and  Marcellus  (1981)  and  of  Wad- 
hams  (1983)  for  a  76-km  line  (the  distance  from 
the  artificial  gravel  island  Kopanoar  to  the  shore). 


The  return  period  for  an  impact  is  taker,  to  be  10, 
100,  and  1000  years.  There  are  large  differences  in 
the  estimates,  with  our  burial  depths  being 
roughly  3  m  less  than  Wadhams.  In  fact,  for  the 
25-m  water  depth  our  estimates  would  only  be 
4.05  and  5.47  m  (assuming  X  =-  3.7  and  2.7  respec¬ 
tively)  if  we  took  g  to  be  20;  a  value  4  times  that 
observed.  We  believe  the  difficulty  with  Wad- 
hams'  approach  lies  not  in  its  principles  but  in  the 
difficulty  o'  obtaining  appropriate  values  to  use  in 
the  thco,'v.  or  instance,  keel  depth  characteristics 
in  deepe:  valet  where  it  is  possible  to  probe  the 
underside  •' r  the  ice  via  submarine  are  probably 
appreciably  -Tiff  ’rent  from  those  in  water  of  50  nr 
or  less  where  gouging  is  currently  taking  place  In 
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addition.  i!  is  at  present  particularly  difficult  to 
know  what  values  to  assume  for  (lie  distance 
dulled  pet  seat  by  the  ice  cos et  user  a  given  point. 
When  gouging  starts,  the  tee  is  slowed  and  mans 
tunes  stopped,  as  the  grounded  ice  tends  to 
stabilize  the  nearby  pack,  converting  n  to  last  ice. 

I  he  ditterenees  between  our  estimates  and  those 
ot  I’llkmgton  and  Mareellus  ( 14S1 )  are  less  by  I  to 
1  m  than  out  differences  witlt  Wadhams’  esti¬ 
mates.  we  tmd  this  somesshat  surprising  as  their 
ptoceduies  appear  to  be  essentially  identical.  The 
ditterenees  m  their  estimates  would  appear  to  be 
largely  the  result  of  ditterenees  in  the  data  used  to 
estimate  the  number  of  gouges  km  yr  .  In  that 
Wadhams  used  direeet  observations  of  keels  while 
I’llkmgton  and  Mareellus  indirectly  inferred  the 
number  ol  keels  from  laser  measurements  of  ridge 
sail',  one  would  expect  Wadhams’  number  to  be 
more  realistic.  C  learly  we  are  a  long  way  from 
achieving  a  concensus  regarding  suitable  pipeline 
burial  depths. 

CONCH  SION 

In  this  paper  we  have  presented  a  large  amount 
of  data  on  the  statistical  characteristics  of  the  ice- 
produced  gouges  that  occur  on  the  Alaskan  shelf 
of  the  Beaufort  Sea  in  shallow  water  (  <  .18  m).  Al¬ 
though  at  first  glance  the  gouges  appear  to  be 
rather  chaotically  distributed,  in  a  statistical  sense 
they  are  very  systematic.  Consequently  we  have 
used  this  information  to  estimate  the  requisite  bur¬ 
ial  depths  of  pipelines  that  would  allow  one  hit  by 
an  ice  mass  in  a  specified  number  of  years. 

In  conclusion  we  would  like  to  comment  on 
some  problems,  the  study  of  which  would  con¬ 
tribute  to  the  understanding  of  the  geophysics  of 
gouging  and  to  the  safe  design  of  sea-floor  pipe¬ 
lines  in  regions  where  gouging  is  known  to  occur. 
We  believe  the  weakest  link  in  the  present  study  is 
the  paucity  of  information  on  the  rate  of  occur¬ 
rence  of  new  gouges  and  their  characteristics. 
Field  programs  should  be  expanded  to  collect  this 
type  of  information.  In  areas  where  offshore 
development  is  contemplated,  it  is  important  to 
start  studies  of  gouging  rates  as  soon  as  possible, 
as  the  collection  of  an  adequate  data  set  takes 
several  years. 

Systematic  regional  sampling  is  also  required  to 
reveal  changes,  if  any,  in  the  probability  density- 
functions  of  parameters  such  as  gouge  depth  with 
changes  in  location  and  in  environment  on  the 
shelf.  Current  information  suggests  that  there  are 


appreciable  changes  in  such  parameter-,  on  a  re 
gional  scale  (lot  instance,  between  the  gouge 
depths  in  the  present  study  area  and  those  ob¬ 
served  otl  the  Mackenzie  Delta).  Studies  should 
also  be  carried  out  to  quantify  the  effects  on  goug¬ 
ing  ot  ditterenees  in  slope  angle  and  aspect  and  the 
nature  of  bed  material.  Such  work,  in  conjunction 
with  detailed  site-specific  studies,  would  be  verv 
useful  m  evaluating  hazards  along  specific  pipeline 
routes. 

Iheoretical  studies  should  also  be  implemented 
to  advance  our  ability  to  treat  gouging  as  a  sto¬ 
chastic  process.  For  instance,  it  would  be  useful  to 
look  at  gouging  as  a  simple  covering  problem  in 
geometric  probability.  If  such  developments  are 
sufficiently  general,  they  can  be  applied  to  differ¬ 
ent  geographic  areas  by  simply  changing  the  val¬ 
ues  of  the  input  parameters. 

We  also  note  that  although  we  have  utilized  an 
exponential  distribution  to  describe  the  relative 
frequency  of  occurrence  of  gouges  of  hffetent 
depths  because  of  its  simplicity  and  the  tact  that 
pressure  ridge  keels  can  be  well  described  by  such 
a  distribution,  \-  tests  of  goodness  of  fit  are  com¬ 
monly  failed.  Therefore  attempts  should  be  made 
to  obtain  a  more  satisfactory  distribution  to  de¬ 
scribe  gouge  depths.  The  same  general  comment 
can  be  made  about  our  utilization  of  a  Poisson 
distribution  to  describe  A  in  that,  as  was  noted 
earlier,  there  are  consistently  more  large  V  10 
values  than  predicted  by  the  fitted  Poisson.  We 
suggest  that  at  least  some  of  these  difficulties  arise 
from  the  fact  that  there  is  no  adequate  treatment 
of  the  infilling  of  the  gouges  in  either  this  or  other 
published  papers  on  gouging.  The  development  of 
a  numerical  simulation  model  that  includes  a  de¬ 
scription  of  both  initial  gouging  and  subsequent 
infilling  of  existing  gouges  could  prove  to  be  illu¬ 
minating. 

Finally,  it  would  be  useful  to  improve  our  un¬ 
derstanding  of  the  interactions  between  pressure 
ridge  and  ice  island  keels  and  the  sea  floor.  Per¬ 
haps  such  studies  will  provide  insight  into  the  pos¬ 
sibility  of  determining  maximum  probable  g<''ige 
depths  for  a  given  sediment  type.  Until  such  mtor- 
mation  is  available,  we  can  only  assume  that  even 
apparently  “impossibly"  deep  gouges  have  a  fi¬ 
nite  probability  of  occurrence. 
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